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A NEW CAPABILITY FOR PREDICTING HELICOPTER ROTOR AND 
PROPELLER NOISE INCLUDING THE EFFECT OF FORWARD 

MOTION 

By F. Farassat* and T. J. Brown 
SUMMARY 


This paper discusses the governing equation and computing technique for 
the prediction of helicopter rotor and propeller noise. The method which 
gives both the acoustic pressure time history and spectrum of the noise includes 
the thickness and the loading noise. It has been effectively adapted to 
computers resulting in a new capability in noise prediction by removing many 
of the restrictions and limitations of previous theories. The capability 
results from the fact that the theory is developed entirely in the time domain 
in contrast to most previous works which were developed in frequency domain. 

The formulation and the technique used is not limited to compact sources, 
steady level flight or to the far-field. In addition, the inputs to the compu- 
ter program are normally available or are amenable to experimental measurements. 
This program can be used to study rotor and propeller uoise with the aim of 
minimizing the radiated noise to reduce annoyance to the public. Several exam- 
ples demonstrating the features and capability of the computer program are 
presented. 


INTRODUCTION 


The problem of noise radiation from propellers and helicopter rotors has 
gained prominence due to its annoyance to the public. Although a helicopter 
rotor is one of the several noise generating sources of helicopters, it is 
most important in the external regions of the present machines. Large propel- 
lers are currently under study for propulsion of large airliners with cruise 
speed of about 850 km/h. The good fuel efficiency of such propulsion systems 
is the main reason for their consideration. Clearly, the reliable prediction 
of the noise of propellers and rotors in the design stage of the aircraft is 
an important step in controlling the level of the noise intensity. 

There has been a steady advance in the last decade in the prediction of 
the noise of rotating blades (ref. 1). There are still disagreements between 
the theoretical and experimental results. The available theories suffer from 
a combination of the following restrictions: 
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a. Compactness of the acoustic sources 

b. Hovering rotor or static propeller 

c. Observer in the far field 

d. Limited airfoil shapes 

e. Limited surface pressure distribution models 

f. Singularities in the solution for high blade tip speeds 

g. Neglect of the thickness noise 

h. Blades with rectangular planforms 

It is believed that the removal of these restrictions and the inclusion of the 
nonlinear propagation effects should result in reliable prediction of the rota- 
ting blade noise. 

Traditionally, rotor noise has been divided into several categories such 
as rotational, vortex and thickness noise. Propeller noise theories consider 
steady thrust and torque forces. However, the acoustic sources are almost 
always assumed compact. Theoretically, the acoustic sources may be grouped 
into two broad classes - those depending on the local pressure and viscous 
stress distribution on the blades and those due to the normal velocity distribu- 
tion on the blades. For example, rotational noise belongs to the first class 
and thickness noise to the second. A theory which incorporates the effects 
of surface pressure and normal velocity distribution on a moving body is devel- 
oped in reference 2. The formulation is then specialized for propellers and 
helicopter rotors. In this work a study of compactness assumption of sources 
on moving bodies has revealed that in the case of helicopter rotors and pro- 
. pellers, the sources on the blades cannot be considered compact for the obser- 
ver position in a large region of space around the rotor. If the compactness 
restriction is removed, then one would like to remove the restrictions of 
limited airfoil > limited airfoil shapes and surface pressure distribution models 
to improve the prediction technique. 

The present paper discusses the governing equation and the computing 
technique together with a computer program developed by the authors at NASA 
Langley Research Center based on the results of reference 2. The purpose of 
developing this program has been to remove the restrictions mentioned above 
and thus improve the prediction of the rotor and propeller noise. The acoustic 
computation is performed in the time domain and the resulting pressure signature 
is then Fourier analyzed to get the acoustic pressure spectrum. 

The new program can only handle deterministic pressure distribution on 
the blades. It is known that the random unsteady pressure fluctuations on the 
blades have significant effect on the noise field of rotating blades when the 
sources can be considered compact. At relatively high tip speeds, particularly 
in forward flight, the unsteady pressure fluctuations are small compared to 
the steady component on propeller blades (ref. 3), At high advancing tip 
speeds, it appears that the random fluctuating pressure of the rotor blades 
has small effect on the acoustic field in the region of space where the acous- 
tic sources should be considered noncompact. It must be mentioned that the 
compactness of sources is not a property of the sources per se but depends 
on the observer position as well as the motion and extent of the sources. 
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Examples are presented in this paper to demonstrate the broad range of 
problems which may be handled by this new program. These examples are selected 
mainly with regard to the restrictions discussed earlier which are removed by 
the new formulation. 
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SYMBOLS 

speed of sound (m/sec) 
constants describing airfoil shape 

blade chord as function of rig (m) 

dimensionless variable n-j/CH 

dimensionless variable ri2/P 

equation describing the surface of each blade 

the equation of a collapsing sphere for a fixed x and t 

integrals given by equations (3) to (7) 

equation of leading edge of the blade as a function of rig 

surface pressure of the blade (N/m ) 

2 

acoustic pressure (N/m ) 

2 

pressure on the lower surface of airfoil (N/m ) 

2 

pressure on the upper surface of tne airfoil (N/m ) 

(Pl Pg)/2 (N/m^) 

|x - yi 

rotor radius (m) 
observer time (sec) 

blade thickness ratio as a function of rio (given as fraction 
equation of trailing edge of the blade as a function of rig 
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a 

r 

Ap 


2 3 

airfoil coordinate given usually as y=T CH(C^E'® + CgE+C^E +C^E 

+05^^) (m) 

source position in a frame fixed to undisturbed medium 
observer position in a frame fixed to undisturbed medium 
blade angle of attack as a function of n 

2 

the curve of intersection of g=0 and f=0 or g=0 and the mean 
plane of each blade 

2 

pressure differential p^^ - py (local lift/unit area) (N/m ) 

coordinate axis fixed to each blade in chordwise direction 
increasing from leading edge to trailing edge 

same as n^-axis but with the origin at the leading edge 
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coordinate axis fixed to each blade in spanwise direction 
with origin at the rotor center increasing towards the 
blade tip 

the angle between radiation direction x - y and the normal 
to blade surface 


pQ density of the undisturbed medium 

T source time (sec) 

T, ,T 2 the source times when the sphere g=0 enters and leaves the 

blade system, respectively (sec) 


THE ACOUSTIC FORMULATION 


The forum! ati on derived in reference 2 is briefly discussed here. Consider 
a moving body whose surface is described by f(^,x)=0 where t is the source time. 
Let v^ be the local normal velocity of the surface, the acoustic pressure 

p'(x,t) is given by 
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The computer program is written in such way that the contributions of the 
upper and lower surfaces of the blades to the acoustic pressure p' are summed 
separately as follows. Let the subscripts U and L stand for upper and lov/er 
surfaces of the blades. If p^ and Ap are defined by the following relations 

|>S “ f <Pl 

Ap = p|^ - Py, (local lift/unit area) 
then, we get 

Pu " Ps " 

Pl = Ps + 


Equation (1) is then written as 
The expressions for I-j to Ig are 
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+ cot 0^_]dr 


dt 


dT 


(7) 


The curve r is now the intersection of the sphere g = r-t+r/c = 0 and the mean 
surface of the blades. The method of computation of the above integrals are 
discussed in the next section. 

The integral describes the thickness noise. The sum of I 2 to will 

be referred to as surface pressure or loa'^ing noise. Aerodynamic calculations 
generally involve the evaluation of Ap. The pressure p^ is usually not 

available. The contribution of the integrals involving p^ to the overall noise 

p'(x,t) appears to be small compared to the dominating term in Eq. (2). 


COMPUTATIONAL METHOD 


Equations (3) to (7) are evaluated on a computer using a double numerical 
integration followed by numerical smoothing and differentiation where required. 
Each of the five integrals are integrated separately. The first three are 
subsequently differentiated and the resulting five pressure contributions are 
added to obtain the pressure signature and spectrum. 

At source t = a sphere is constructed with its center at the observer 

location. Its radius r-j is selected such that its circle of intersection, 

C' , with the plane of the rotor is tangent to the rotor disk. From this 
initial geometry the initial observer time, t^ , is calculated from t^. = x^. + 

r^c where c is the speed of sound in the medium. The sphere is allowed to 

collapse by an amount cAx, where x is the emission or source time. During 
this period, the helicopter rotor is allowed to translate and rotate. The 
resulting arc of intersection between the rotor disk and the new C is swept 
point by point in a counterclockwise direction until an intersection with a 
blade surface is detected or until the arc passes out of the rotor disk. When 
a blade is encountered, the integrands of equations (3) to (7) are evaluated 
and subsequently the line integrals are accumulated point by point using a 
trapezoidal scheme. 

The collapsing process of the sphere g = 0 is repeated, each time yeilding 
a value for the line integrals which are accumulated for the source time 
integration using simpson rule. This process is continued until it is detected 
that the collapsing sphere has passed out of the rotor disk. The integration 
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is thus concluded for the oberserver time t. and the resulting integrals are 
saved for further processing. Successive points are obtained in like manner. 

To facilitate numerical smoothing and differentiation with respect to the 
observer time t, it is required that the t^'s be equally spaced. Since the 

relation between the oberserver time t and the source time t is in general 
nonlinear, an iteration technique is used to obtain the initial radius r^ and 

the corresponding source time where the sphere g = 0 begins to collapse. The 

smoothing and numerical differentiation which is used are presented in reference 
4. It is based on the theory of finite Fourier series using factors which 
modify harmonic levels to improve convergence characteristics and to reduce 
Gibbs phenomenon. As a byproduct of this, the pressure spectrum of the 
acoustic signature is obtained quite easily using intermediate results of the 
smoothing and differentiation process. 

The observer may be assumed fixed in the moving frame attached to the 
vehicle. In this case two-point time differentiation is used while the observer 
position is frozen with respect to the undisturbed medium. For a new observer 
time, the observer is moved to its original position in the moving frame and 
the process is repeated. 


SOME APPLICATIONS OF THE COMPUTER PROGRAM 


The examples in this section are selected with realistic data to demonstrate 
the features and capability of the computer program. In all the calculations, 
the source distribution is noncompact. In each example, the assumptions and 
the needed data used in computation are given. In two cases, experimental 
measurements are also available and are presented together with theoretical 
results. 

EXAMPLE 1 - Helicopter noise 

The special features of this example are; 
i) helicopter in flight, 

if) triangular blade tips and blades with twist, 
iii) realistic rotor attitude - tip path plane does not contain the 
forward velocity vector. 

Because of the high advancing blade tip speed and the observer location, only 
thickness noise is believed to be dominant and therefore is calculated. The 
pressure distribution on the blades was not available. 
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INPUT DATA: 


Number of blades = 2 
R = 7.62 m 
RPM = 300 


CH(r)2) = 


0.838 m 

-2.246 ri2+17.283 


T(n2) = 
LE(n2) = 


0.1 

-0.167802 + 1.3286 
-0,335 m 

. 2.24602 - 16.780 m 


025.7.322 m 
m 7. 322 <0217. 62 m 


025.7.322 m 
7.322<02l7.62 m 

0 2 1 7.322 m 
7.322<02<7.62 m 


TE(o 2) = 0.503 m 

a(o 2 ) = -8.0+1.0502 degrees 

y(o]>n 2 ) = T(o2)CH(o 2)(3.3333 E - 6.5079 E^ + 3.1746 E^} m 


c = 340. m/sec 

Helicopter speed = 259.3 km/h (140 kots), level flight 
Helicopter altitude = 152.5 m (500 ft) 

Observer positioo (at emissioo time) = 439.6 m (1441.3 ft) directly ahead 
of helicopter at grouod level (Fig. 1). 

Rotor angle of attack = 6.5 degrees 

Figure (1) shows the calculated and experimental acoustic pressure signatures. 
The experimental pressure signature includes the tail rotor noise. The 
helicopter speed and altitude used in this example are those recorded from the 
aircraft instruments. The emission distance is, however, an estimate obtained 
from approximate overhead position of the helicopter. The test helicopter has 
blades with Wortmann FX69-H-098 airfoil section which has a small leading edge 
radius. The airfoil section used in calculations is biconvex cubic shape 
approximating the actual airfoil section. Based on manufacturer supplied noise 
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data, the calculated sound pressure level is good. As seen from Figure 1, the 
agreement between the shape of the theoretical and measured signature is also 
good. 

EXAMPLE 2 - Propeller noise 

The special features of this example are: 

i) blades with twist and variable thickness 

ii) observer in motion with the aircraft 

The power absorbed by the propeller was available in this case and the pressure 
distribution Ap was used in the calculation. This distribution is thought 
to be realistic. To simplify the calculation of this distribution, the blade 
planform is assumed rectangular. The actual blades on the test aircraft had 
a slight taper. The blade form curves are presented in Fig. 2. The thickness 
ratio and blade twist are approximated using the curves in this figure. 

INPUT DATA: 

Number of blades = 3 

R = 1 . 30 m 

RPM = 2145 

CH = 0.156 m (uniform, taken as chord at .85 R) 

TCng) = 0.069 + 3.2244 exp(-8.615 ^^) 
airtg) = 3.61 + 78.037 exp (-2. 685 rig) degrees 
y(n,',%) = CH T(nn) (3,3333 E - 6.5079 E^ + 3.1746 E^) 
c = 345 m/sec 

Aircraft speed = 144.5 km/h (78 KTS) 

Observer position (in disk plane) = 7.28 tn from propeller center moving with 
aircraft. 

Ap(n^,ri 2 ) = 2.507 x 10^ E*^^^(l-E*^) E^-®^^(l-E^) N/m^ 

Figure (3) shows the calculated and measured acoustic pressure signature 
and spectrum. The agreement both in time domain and frequency domain is very 
good. Note that although steady source distribution is used, the high har- 
monics of the acoustic pressure spectrum are calculated with reasonable 
accuracy. For propellers in motion, the blade geometry and steady blade 
surface pressure distribution are believed to be sufficient for prediction of 
the generated noise. 


9 



EXAMPLE 3 - Rotating Blades 


The special features of this example are: 

i) supersonic tip speed 

ii) The blade surface pressure distributions Ap and p^ are included 
in noise calculation 

To study the contribution of the terms involving p^ to the acoustic pressure, 

this example was devised. In general, p^ is not available and is not calculated 
usually. This is because p^ does not contribute to the lift, thrust, or the 

torque on the engine. Linearized aerodynamic theory is used to calculate Ap 
and p^ assuming that at each radial position, the flow around the blade is 

two dimensional. The flow around the entire blade is supersonic. A uniform 
downwash of 40 m/sec with tip Mach number of 1.375 are assumed. The blade 
Mach number at the inner radius is 1.10. The blade disk is assumed to be 
static. 


INPUT DATA; 


Number of blades = 2 
R = 5.0 m 

Inner blade radius = 4.0 m (Blade length = Im) 


RPM = 906.0 
CH(ri 2 ) = 0.4 m 
T(h2) = 0-05 

a(r,2) = 


0.768v^Rpr 


degrees, 


\ " 02 ^ 2 /^ (Local Mach Number), 

= angular velocity of the blades 
(0.05n' m 0£n'£-2 m 

y(n' ) =S ] ^ 

1 2 (0.02 (l-2.5n^) m 0.2<n^£.4 m (diamond section) 

Observer position (in plane) = 15 m from disk center 

Observer position (out of plane) = 30° above disk plane and 15 m from disk 
center. 
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Ap(n',n ) = 3840.0 N/m^ 
1 2 

' 7.159 X 1Q3 
-7.159 X 10^ 


Ps(n^.n )= 


N/tn2 


0 <n^< 0,2 m 

0.2<n'<0.4 m 
1 


Figure (4) presents the separate contributions of each term in Eq. (2) 
to p'('^,t} together with p'(x',t) both in time and^frequency domains. For 
both observer positions, the acoustic pressure p'{x,t) is dominated by the 
thickness noise. The shapes of the acoustic pressure signature and spectrum 
for the observer in the disk plane are substantially different from those for 
the observer above the plane. This signifies a complex wave structure 
rotating in or near the plane containing the blade disk (not necessarily in 
the near field only). 


In the plane of the disk, the contribution of the term involving Pg is 

larger than that due to Ap. These contributions are of the same phase. 

Above the plane this trend reverses although the contribution of Pg term is 

still relatively large. Since the thickness noise is significantly higher 
than other contributions, the inclusion of terms Ig to Ig produces small 

changes in the acoustic pressure signature and spectrum. However, for high tip 
speeds, if the effect of Ap is included in the calculations, one should also 
include the effect of Pg. Some calculations at subsonic tip speeds have shown 

that the contributions of the terms involving P 5 is smaller than the other 

terms. However, this requires further study. 

EXAMPLE 4 - (Helicopter rotor) 


The special feature of this example is blades with swept back tips. Only 
thickness noise is calculated. This example is worked out to show the favor- 
able effect of nonrectangular blade planform. 


INPUT DATA 


Number of blades = 2 


R = 5.. m 
RPM = 527.1 

CH(ri 2 ) = 0*4 m (see Fig. 5) 

T(ri2) " 

(- 0.2 m x\^<(\.2.S m 

'"^^'^2 ' 1-3.019 + 0.6633 r\^ m 4.25<n2l5 m 
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TE{ri2) = 


0.2 m ri2<4.25 m 

•*2.619 + 0.6633 TI 2 m 4.25<n2*^5 m 


ot(ri2) = 0 

y(n^ , 712 ) " 0.064 (E-E^) biconvfix parabolic 
c = 346 m/sec 

Helicopter speed = 248.4 km/h (134.1 kts) level flight 

Observer position (in rotor plane) = 50. m from rotor center at the start of 
emission 


Figure 5 presents the calculated acoustic pressure signature and spectrum 
(thickness noise). For comparison the corresponding results for a rotor with 
rectangular blades (CH = 0.4 m) are presented. For the rotor with rectangular 
blades, all the input data (except, of course, for TE(n 2 ) LE(ri 2 )) 

identical to those given above. It is seen that the sweep at the blade tip 
has considerable influence in the shape and the peak values of the acoustic 
pressure signature. As compared to that of rectangular blades the acoustic 
pressure spectrum shows reduction in the case of swept back blades up to 
the 30th harmonic. The maximum reduction is about 6 dB at about 15th harmonic. 
This example demonstrates that the planform variation should be considered a 
promising method of controling the noise of helicopter rotors and propellers. 


CONCLUSIONS 


This paper presents a new formulation and a discussion of the method of 
computation of helicopter rotor and propeller noise. Only deterministic 
time dependent blade surface pressure may be used in the computer program that 
has been developed based on the new fonnulation. There are many situations 
where the unsteady random pressure fluctuations do not contribute substantially 
to the acoustic pressure. The most common of these is in the case of rotating 
blades at high tip speeds. The examples in this paper demonstrate the range 
of applicability of the computer program. By removing the restrictions and 
limitations of previous theories, it provides a capability which will improve 
the prediction and reduction of rotor and propeller noise. 

Random unsteady pressure fluctuations are important in regions where the 
acoustic sources are compact. They appear to be also important for static 
propellers and hovering rotors even in the regions where the acoustic sources 
are noncompact [ref. 3]. This is because of the injestion of atmospheric 
turbulent eddied. In this case, a combined noncompact source calculation 
for thickness and steady loading noise and compact source calculation for 
unsteady loading noise appears to be the best choice. 
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(a) Helicopter position at (b) Blade Tip Geometry 

emission time 
HO = 439.6 m 
TPP: Tip Path Plane 



Figure 1 (Example 1) - Emission Distance, Blade Tip Geometry and Comparison 
Between Theory and Experiment for a Helicopter in Flight. 
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